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ABSTRACT 

We compare the improved measurement of the Lya forest flux probabihty distribution 
at 1.7 < z < 3.2 presented by Kim et al. (2007) to a large set of hydrodynamical 
simulations of the Lya forest with different cosmological parameters and thermal 
histories. The simulations are in good agreement with the observational data if the 
temperature-density relation for the low density intergalactic medium (IGM), T = 
TqA^"-'^, is either close to isothermal or inverted (7 < 1). Our results suggest that the 
voids in the IGM may be significantly hotter and the thermal state of the low density 
IGM may be substantially more complex than is usually assumed at these redshifts. 
We discuss radiative transfer effects which alter the spectral shape of ionising radiation 
during the epoch of He 11 reionisation as a possible physical mechanism for achieving 
an inverted temperature-density relation at z ~ 3. 
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1 INTRODUCTION 

Traditional analyses of the Lya forest observed in the spec- 
tra of high redshift quasi-stellar objects (QSOs) decom- 
pose the flux distribution into a series of discrete absorp- 
tion profl l es, Renerally characte r ised by a Voigt function 



iIHu et al 
iKim et al 



iLidz et ai] I2OO6I: [Becker et all 120071') and wavelet analysis 
l|MeiksiDll2000l : iTheuns fc Zaroubil I2OO0I : IZaldarriagal I2OO2I ) 
have also been developed. These measures are simple and - 
at least in principle - easy to compare to the same quan- 
titie s extracte d from theoretical models of the Lya for- 
est (|Cenlll992l: iMiralda-Escude et allll996l: IbI fc DayidsenI 



I995I : iLu et al.l 1 19961 : iKirkman fc Tvtled 1 19971 : | l997l:ICroft et al.lll999l : lDave et al.ll 19991 : iTheuns et al.lll998l : 



19971 . [2OOII). Voigt profiles provide an accurate 



description for an absorption line if the absorber i s a lo- 
calise d gas cloud with a Gaussian velocity dispersion (jRauchl 
119981 ). In the now widely established paradigm for the 
origin of the Lya forest most of the absorption is, how- 
ever, caused by extended gas distributions broadened by 
the Hubble flow, removing much of the physical motiva- 
tion for the decomposition of the flux distribution into 
Voigt profiles. Such a decomposition is nevertheless use- 
ful, but it is non-unique and the line fitting process is 
time consuming and somewhat subjective (see iKim et al.l 
I2OO7I . hereafter K07, for an exhaustive discussion). Al- 
ternative characteri sations of the flux distr i bution based 
on p ix el statistics ([Jenkins fc Ostrikerl 199ll: Ranch et al.l 
19971: iGaztanaga fc CrofiJ Il999l: [ McDonald et all 1200' " 



Jena et al.ll2005h . In particular, the power spectrum of the 
Lya forest flux distribution has recently been successfully 
developed into a quantitative tool for measuring the matter 
power spectrum on scales of 0.5fe~^ Mpc to AQh~^ Mpc (e. g. 
ICroft et al.|[2002l : I Viel et al.|[2004bl : [McDonald erallboOSl ). 

The simplest pixel statistic is the Lya flux probability 
distribution function (PDF) , which is sensitive to the density 
distri bution and thermal state of the IGM (e.g. iBecker et al.l 
[2007]). Several attempts have been made to obtain joint con- 
straints on a variety of cosmological and astrophysical pa- 
rameters using the flux PDF togeth e r with the flux power 



spectrum (ICh oudhurv et al.l 2001 ; Desiacques fc Nusseil 



Theuns et all I2OO0I : iMeiksin et al.l I2OOII : IViel et al.l 12004^ 



I2OO5I : [Lidz et al..,2006, : .Desiacques et al.|[2007l '). However, the 
PDF is also sensitive to a range of systematic uncertainties, 
most notably the assumed continuum level, the noise prop- 
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erties of the spectra and contamination by metal absorption 
lines in the Lya forest. In high resolution, high signal to 
noise data at 2; < 3, the first two uncertainties have a small 
impact on the shape of the PDF (K07). The 2 < 3 metal ab- 
sorption in the Lya forest, however, alters the shape and am- 
plitude of the PDF significantly, especially at 0.2 < F < 0.8. 
Since many numerical simulations of the Lya forest do not 
include metals, and those that do must also correctly repro- 
duce their uncertain distribution and abundance, this can 
lead to substantial systematic uncertainties when attempt- 
ing to fit simulation data to the observed flux PDF of the 
Lya forest. 

K07 have completed a detailed analysis of the impact 
of metal absorption lines on the Lya PDF at 1.7 < z < 3.2. 
They presented new measurements of the PDF in which 
all identifiable metal contamination has been removed from 
the Lya forest on a line-by-line basis. The metal cleaned 
PDF had a significantly different shape when compared di- 
rectly to the metal contaminated sample. In this paper we 
shall investigate what can be inferred about the underly- 
ing matter distribution and the thermal state of the inter- 
galactic medium from these new PDF measurements using 
state-of-the-art hydrodynamical simulations of the Lya for- 
est. Our approach is different to most recent studies of the 
flux PDF which have used either pseudo-hydrodynamical 
techniq ues, dark matter only simulations or sem i-analytical 
models ( Choudhurv et al.l 2001 : Desiacaues fc Nusser 2005. 



iLidz et al.ll2006l : lBecker et al.ll2007l ') which model the effect 
of gas pressure on the underlying IGM density distribution 
in an approximate way. 

The structure of this paper is as follows. We discuss the 
observational data and the hydrodynamical simulations in 
section 2. In section 3, we qualitatively compare the simu- 
lation data to the observations and determine the impact 
of various astrophysical and numerical effects on the sim- 
ulated flux PDF. A quantitative comparison is undertaken 
in section 4 and the best fitting models to the K07 data 
are presented along with a comparison to other recent stud- 
ies in the literature. We briefly discuss the implications of 
our results for measurements of the matter power spectrum 
from the Lya forest flux power spectrum in section 5 and 
conclude in section 6. 



2 OBSERVATIONS AND SIMULATIONS OF 
THE Lya FOREST 

2.1 The observational data 

The Lya flux PDF at 1.7 < z < 3.2 was recently measured 
by K07 using a set of 18 high resolution {R ~ 45 000), 
high signal to noise (S/N > 30 - 50) VLT/UVES spectra. 
Metal lines were removed from the Lya forest by directly 
fltting all identifiable absorption features. This approach to 
metal absorption removal is different to other methods in 
the literature, w hich either account for metal contami na- 
tion statistically (iTvtler et al.ll2004lK[rkman et al.ll2005l ) or 
by excising heavily contaminated regions from the spectr a 
JRauch et al.lll997l : iMcDonald et al.ll2000l : iLidz et al.ll200^ . 
Following this procedure, K07 obtained the Lya flux PDF in 
three redshift bins at (z) = 2.07, (z) = 2.52 and {z) = 2.94. 
We shall compare these new improved measurements to the 



PDF computed by numerical simulations of the Lya forest. 
Further details regarding the observational data and its re- 
duction may be found in K07. 



2.2 Numerical simulations of the IGM 

The simulations used in this stud y were run usin g the par- 
allel Tree-SPH code gadget-2 (|Springell |2005^ . modified 
to compute the ionisation state of the gas using the non- 
equilibrium ionisation algorithm of Bolton et al. (2006). The 
fiducial simulation volume is a periodic box 20/i^^ comoving 
Mpc in length containing 2 x 256^ gas and dark matter par- 
ticles. The gravitational softening length was set to 1/30"' 
of the mean linear interparticle spacing. The gas in the sim- 
ulations is photoionised by a spatially uniform ultraviolet 
background ( UVB) based on the gala jcies and QSOs emis- 
sion model of lHaardt fc Madaru (|200H ). The mean free path 
for hydrogen ionising photons is expected to be much larger 
than the mean separation between ionising sources at 2 < 3, 
and so the assumption of a uniform U VB at the hydrogen 
Lyman limit should be a reasonable one (jBolton fc Haehneld 
l2007bl ). Star formation is included using a simplified pre- 
scription which converts all gas particles with overdensity 
A > 10'' and temperature T < 10^ K into coUisionless stars, 
significantly speeding up the simulations. The runs were all 
started at 2 = 99, with initial condition s gene rated using 
the transfer function of lEisenstein fc Hul l|l999l ). The sim- 
ulations explore a range of cosmological and astrophysical 
parameters consistent with other observational data, and are 
summarised fully in Table 1. 

Constraints on the thermal state of the low density IGM 
are usually expressed in terms of a polytropic temperature- 
density relatioiu, T = TqA^^^, expected to arise for A — 
p/{p) ^ 10 when photo- heating by t he UVB is balanced 
by cooling due to ad iabatic expansion (JHui fc Gnedinlll997l : 
IValageas et al.l I2OO2I ). Measurements of the thermal state 
of the IGM at 1.7 < 2 < 3.2 from Lya forest data 
give temperatures in the range 10^'^ K < To < 10*'^ K 
and are consistent with slopes in the range of those ex- 
pected within the s tandard paradigm f or the photoheated 
IGM, 1 < 7 < 1.6 jRicotti et al.l I2OO0I : ISchave et aj] I2OO0I : 
iMcDonald et al.ll200ll '). The fiducial model (20-256) was cho- 
sen to be representative of these constraints. In order to ex- 
plore the effect of different temperature-density relations on 
the flux PDF, we also run six further simulations with differ- 
ent thermal histories. Four simulations (20-25611, 20-25612, 
20-256gl and 20-256g2) reproduce values of To and 7 consis- 
tent with the upper and lower end of the range quoted above. 
We have also simulated two further models with 7 < 1 
(20-256g3) and an additional scatter in the temperature- 
density relation (20-256xs). We shall discuss the motivation 
behind these latter choices further in section 3. The dif- 
ferent thermal histories are constructed by modifying the 
fiducial simulation He 11 photo-heating rate, eHoiii such that 
ehcII = aA''ejj'Jii. In the case of the model with extra scat- 
ter (20-256xs), a random Gaussian distributed dispersion 



^ The temperature-density relation is also sometimes referred to 
as the IGM effective equation of state, or simply just the equation 
of state. Note, however, that the equation of state for this low 
density gas in the correct sense is the ideal gas equation. 
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Table 1. Parameters used in our suite of hydrodynamical simulations. In all models we assume a flat Universe with Ha = 1 — f!m, 
f!m = 0.26, f2b/i = 0.024, h = 0.72 and n = 0.95. The temperature-density relation is estimated at each redshift by a least squares fit to 
the volume weighted temperature-density plane at O.f < A < 1 and To < 10^ K. The To values have been rounded to three significant 
figures. 



Model 


o"8 


Tom 

(z = 2.07) 


To[K] 
(z = 2.52) 


To[K] 
(z = 2.94) 


7 
{z = 2.07) 


7 
(2 = 2.52) 


7 
{z = 2.94) 


a 


/9 


Notes 


20-256 


0.85 


14800 


17600 


20800 


1.37 


1.31 


1.27 


1.0 


0.0 


Fiducial model 


20-fOO 


0.85 


16800 


19600 


22700 


1.40 


1.35 


1.30 


1.0 


0.0 


Low resolution 


20-200 


0.85 


15100 


18000 


21300 


1.38 


1.33 


1.29 


1.0 


0.0 


Med. resolution 


20-400 


0.85 


14400 


17200 


20300 


1.37 


1.31 


1.26 


1.0 


0.0 


High resolution. 


40-200 


0.85 


16900 


19700 


22700 


1.42 


1.36 


1.31 


1.0 


0.0 


Low res. /Med. box 


80-400 


0.85 


16900 


19800 


22800 


1.41 


1.36 


1.31 


1.0 


0.0 


Low res. /Large box 


20-256m 


0.85 


14800 


17700 


20900 


1.37 


1.31 


1.27 


1.0 


0.0 


Multiphase SF 


20-256w 


0.85 


18900 


20900 


23300 


1.47 


1.39 


1.32 


1.0 


0.0 


Strong winds 


20-256tl 


0.85 


10900 


12600 


14600 


1.41 


1.36 


1.31 


0.5 


0.0 


Low To 


20-256t2 


0.85 


17600 


21300 


25500 


1.35 


1.29 


1.25 


1.4 


0.0 


High To 


20-256gl 


0.85 


15300 


18300 


21700 


1.59 


1.56 


1.53 


1.0 


0.5 


High 7 


20-256g2 


0.85 


14100 


16700 


19800 


1.16 


1.08 


1.02 


1.0 


-0.4 


Low 7 


20-256g3 


0.85 


12600 


14700 


17200 


0.67 


0.54 


0.44 


1.0 


-1.2 


Inverted 7 


20-256XS 


0.85 


15900 


18800 


21800 


1.38 


1.32 


1.27 


1.0 


0.0 


Extra scatter 


20-256sf 


0.75 


14600 


17400 


20700 


1.37 


1.32 


1.27 


1.0 


0.0 


0-8 = 0.75 


20-256s2 


0.95 


15000 


17700 


21000 


1.37 


1.31 


1.26 


1.0 


0.0 


erg = 0.95 


20-256s2g2 


0.95 


14200 


16700 


19700 


1.15 


1.07 


1.01 


1.0 


-0.4 


as = 0.95, Low 7 



is applied to the He 11 photo-heating rate. In this way we 
generate self-consistent thermal histories, in the sense that 
the thermal state of the IGM is correctly coupled to the 
gas pressure. Note, however, that a fully physically consis- 
tent treatment of the IGM thermal history at 2 ~ 3 would 
require cosmological radiative transfer simulations with a 
level of detail beyon d at least our current numerical capa- 
biliti e s (although see Bolton et al.l |2004| : iMaseUi fc Ferraral 
I2OO5I : iTittlev fc Meiksinll2007l ). The values of a and f3 used 
in all the simulations are summarised in Table 1, along with 
the resulting parameters for the temperature-density rela- 
tion in each simulation, [j 

The cosmological parameters used for our fiducial sim- 
ulation (20-256) are consistent with the combined anal- 
ysis of the third year WMAP and Lya forest data. 
(nIn,^A,^b/^^fe,c^8,n) = (0 26,0.7 4, 0.024, 0.72, 0.85, 0.95) 
jViel et al.ll2006l : ISeliak et al.ll2006bl ). We also run three ad- 
ditional simulations with different values for ag : one which is 
consistent with the third year WMAP data alone (20-256sl, 
[SpcrRcl ot al. 2007) and another two consistent with the 
somewhat higher constraints from weak lensing and LyQ for - 
est data (20-256s2 and 20-256s2g2. lLesgourgues et al.ll2007l ). 
The third model (20-256s2g2) also has a slope for the 
temperature-density relation which is consistent with the 



^ The values of 7 in the fiducial model are smaller than the 
value of 7 ~ 1.6 typically reached under the assumption that 
the low density gas is optically thin and in photoionisation 
equilibrium following reionisation. We have not rescaled the 
He II photo-heating rate in this instance. This difference is due 
to non-equilibrium ionisation effe cts which occur during rapid 
changes in the UV background (e.g. J Mi ralda- Escude fc Reesll994l : 
iHaehnelt fc SteinmetJ 19981: iTheuns et al.ll998^ . The smaller val- 
ues of 7 are due to the hardening of the UV background spectrum 
around z ~ 3.5 due to the increased contribution of QSOs to the 
metagalactic UV background, leading to the reionisation of the 
He II in the simulation. 



Table 2. Mass resolution and box size of the five additional simu- 
lations used to check numerical convergence. The mass resolution 
and box size of the fiducial simulation is also listed for comparison 
at the bottom of the table. 



Model 



Box size 
[comoving Mpc/h] 



Total particle 
number 



Gas particle 
mass [Mq/H] 



20-100 


20 


2 X 100^ 


1.03 X 10* 


20-200 


20 


2 X 200^ 


1.29 X 10'' 


20-400 


20 


2 X 400=^ 


1.61 X 10^ 


40-200 


40 


2 X 200=^ 


1.03 X 10* 


80-400 


80 


2 X 400^ 


1.03 X 10* 



20-256 



20 



2 X 256=* 



6.13 X 10*5 



lower end of the observational constraints. Two further mod- 
els were_jiHi_jasing_th«niult^ star formation model 
of ISpringel fc Hernguistl (l2003al ). One model includes the 
multi-phase star formation with winds disabled (20-256m), 
while the other includes the effect of galactic winds with a 
velocity of 484 km s^^ (20-256w). This latter model is ex- 
treme in the sense that the kinetic energy of the winds is 
compar able to the energy output fr om supernovae in the 
model (ISpringel fc Hernguistl l2003bl ) and should therefore 
provide an upper limit on the impact galactic winds may 
have on the PDF. 

Lastly, to check numerical convergence we run five ad- 
ditional simulations with different box sizes (40-200, 80-400) 
and mass resolutions (20-100, 20-200, 20-400). The param- 
eters of these simulations are listed in Table 1 along with a 
summary of all the other models, while the parameters char- 
acterising their resolution and box-size are listed in Table 2. 
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2.3 Synthetic Lya spectra 

Sets of synthetic Lya spectra are constructed a.t z — 2.07, 
z — 2.52 and z — 2.94 using 1 024 random lines of sight 
drawn from the simulations (e.g. iTheuns et al.lll998r ). The 
line profile convo lution is performed w ith the Voigt profile 
approximation of lXepper Garcial l|200q '). We rescale the syn- 
thetic spectra to match the mean normalised flux, (F) , in the 
redshift bins used by K07 to measure the flux PDF; (_F) = 
[0.863,0.797,0.730] at (z) = [2.07,2.52,2.94]. These values 
correspond to an effective optical depth, Toff — —\n{{F)) = 
[0.147,0.227,0.315]. 

The raw spectra are then processed to resemble the ob- 
servational data. The spectra are convolved with a Gaussian 
with FWHM = 7 km s~^ and are rebinned onto pixels of 
width 0.05 A. Noise is then ad ded t o the spectra in a simila r 
manner to lRauch et"aLl (|l997l ') and lMcDonald et al.l l|2000h . 
Within each of the three observationally deflned redshift 
bins for the PDF, the variance of the noise in the observed 
spectra is determined over flux ranges corresponding to each 
flux bin of the PDF. Not e that we have used the same bin- 
ning for the flux PDF as iMcDonald et all (|2000l ') and K07. 
The PDF bins have a width of AF = 0.05, with the first 
and last bins centred on F = and F = 1 respectively. This 
gives 21 bins in total from F = to F = 1. Gaussian dis- 
tributed noise consistent with the variances in each flux bin 
of the observational data is then added to the corresponding 
pixels in the synthetic spectra. This process mimics the noise 
properties of the observed spectra well (Ranch et al. 1997). 
Note that all pixels with flux levels smaller than F = 0.025 
or greater than F — 0.975 have been allocated to the F = 
or F = 1 PDF bins, respectively. 



3 COMPARISON OF THE K07 PDF TO 
HYDRODYNAMICAL SIMULATIONS 

3.1 The effect of resolution, box-size and galactic 
winds 

Fig. [l] displays the effect of resolution, box-size and galactic 
winds on the simulated flux PDF. The observational data of 
K07 are shown by the open diamonds with error bars, and 
the shaded regions in all panels correspond to the flux range 
over which noise properties and continuum placement have 
a smaU impact on the PDF, 0.1 < F < 0.8 (K07). 

The top row in Fig. [l] displays the effect of mass resolu- 
tion on the simulated PDF. The ratio of the three simulated 
PDFs with different mass resolutions to the fiducial model 
PDF (20-256, see Table 2 for details) is shown in the lower 
section in each of the panels. The simulated PDF appears 
to be marginally converged and is generally to within a few 
percent of the highest resolution run at all redshifts. As we 
shall see in the next section, these differences are small com- 
pared to those exhibited when other parameters are varied, 
most notably erg and 7. Note, however, the flux PDF of the 
fiducial model does not agree well with the PDF from the 
observational data of K07. We shall return to this point later 
in this paper. 

In the middle row, the effect of the simulation box size 
on the PDF is displayed. The lower sections of each panel 
show the ratio of the PDFs from the two smaller boxes to 
the PDF extracted from the largest box (80-400). All three 



simulations have the same mass resolution. Again, the dif- 
ferences in the simulated PDFs are very small. We therefore 
consider our fiducial simulation box size and mass resolution 
to be adequate for distinguishing between the various IGM 
models considered in this work. 

The bottom row of panels in Fig. [T] shows a comparison 
of the PDFs computed from the fiducial model (20-256, solid 
curve) to the results from the simulati ons including the mul- 
tiphase star formation (SF) model of ISpringel fc Hernguistl 
(|2003al ') with (20-256w, dotted curve) and without (20-256m, 
dashed curve) winds. The lower sections in each panel again 
display the ratio of these PDFs to the fiducial model. The 
differences between the fiducial and multiphase SF model 
PDFs are within a few percent, a result which gives us con- 
fidence the simple star formation prescription used for the 
majority of the simulations in this work has little impact 
on the properties of the Lya forest. This is because the 
highly overdense regions containing most of the star for- 
mation correspond to only a very small proportion of the 
volume probed by the synthetic spectra. The difference be- 
tween the model with strong galactic winds is also small, ex- 
cept at _F = in the (z) — 2.07 bin. However, this is a rather 
extreme model and may overestimate the impact of winds on 
the Lya forest. In compari son, after resc a ling th e mean ffux 
of their synthetic spectra, iTheuns et al.l (|2002cr ) find winds 
also have a very small effect on th e Lya f orest. A similar re- 
sult w as found by iBruscoli et al.l l|2003l ') . iBertone fc Whitg 
(J2006I ) used a semi-analytical model to investigate the im- 
pact of galactic winds on the Lya forest, and also found the 
impact of winds on the PDF is negligible. Our results are in 
qualitative agreement with these studies, and we expect the 
effect of galactic winds on the Lya flux PDF to be small. Al- 
though we have not analysed the possible impact of outflows 
from rarer, more massive galaxies hosting active galactic nu- 
clei (AGN) on the Lya forest, these are expected to extend 
into a smaller v olume of the I GM than the supernova driven 
galactic winds (|Siiacki et al.ll2007f ). AGN feedback is there- 
fore also likely to have a negligible impact on the Lya forest 
flux distribution. 



3.2 The effect of cosmological and astrophysical 
parameters 

The flux PDFs constructed from the simulations with differ- 
ent cosmological and astrophysical parameters are compared 
to the observational data of K07 (open diamonds with er- 
ror bars) in Fig. [2] From left to right, in each of the four 
panels the effect of varying To, as, 7 and Teff on the PDF 
is shown. The solid curve in each panel corresponds to the 
PDF of the flducial model (20-256). The lower segment of 
each panel shows the ratio of the PDFs indicated in the di- 
agram to the fiducial model. Note again that except where 
indicated in the rightmost panels of Fig. (2] all the synthetic 
spectra are rescaled to have t^s = [0.147,0.227,0.315] at 
(z) = [2.07, 2.52, 2.94], corresponding to the effective optical 
depth in the redshift bins for which K07 have reported their 
measurement of the flux PDF. 

At all redshifts the PDF is relatively insensitive 
to changing the IGM temperature at mean density. To, 
for values within the observationally constrained range 
(jSchave et al.l 120001 ). Higher temperatures slightly increase 
the amplitude of the PDF at 0.1 < F < 0.9 and decrease 
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Figure 1. The impact of resolution, box-size and galactic winds on the simulated flux PDF. Top row: The effect of mass resolution on 
the simulated flux PDF. The lower section in each panel displays the ratio of the flux PDFs to the fiducial model (20-256). Middle row: 
The effect of box size. The lower section in each panel displays the ratio of the flux PDFs to the model with the largest box size (80-400). 
Bottom row: The effect of the adopted star formation model and galactic winds. The lower section in each panel again displays the ratio 
of the flux PDFs to the flducial model (20-256). For comparison, the observational data of K07 are also shown (open diamonds with error 
bars) and the shaded grey region in each panel corresponds to the portion of the PDF which is least affected by noise and continuum 
uncertainties. 
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Figure 2. Comparison of tlic flux PDF computed from liydrodynamical Lyo forest simulations witli different cosmological and astropfiysical 
parameters to the observational data of K07. Top row: The observed PDF at [z) = 2.07 (open diamonds with error bars). The solid curve 
in each panel shows the PDF derived from the flducial model. From left to right, the dotted and dashed curves demonstrate the effect on 
the PDF of varying To, erg, 7 and Toff within the observationally determined boundaries. The lower section of each plot shows the ratio 
of these PDFs to the fiducial model. Middle Row: As for top row but at (2) = 2.52. Bottom Row: As for top row but at {z) = 2.94. The 
shaded grey region in each panel corresponds to the portion of the PDF which is least affected by noise and continuum uncertainties. 



© 0000 RAS, MNRAS 000, 000-000 



The flux distribution of the Lya forest 7 



the amplitude at F = and F — 1. This behaviour is 
attributable to the Jeans smoothing effect of gas pressure 
on the spectra; higher temperatures increase the IGM gas 
pressure and smooth the IGM over larger scales. Combined 
with the increased Doppler widths in the higher temperature 
models, this redistributes the PDF by effectively erasing the 
small scale power in the spectra. Varying as has a slightly 
larger effect. Models with smaller values of ag tend to pro- 
duce broader absorption lines which arise from the slightly 
smoother gas distribution at scales comparable to the Jeans 
length. The effect is to again erase small scale power. The 
effect of erg and temperature on the flux PDF are, however, 
small and not sufficient to account for the differences be- 
tween the K07 data and the simulations. A more detailed 
discussion of th e effect of these param eters on flux statistics 
can be found in lTheuns et al.l l|200d ). 

Varying the temperature-density relation index, 7, has 
a much stronger effect on the PDF. Altering 7 changes the 
temperature of the IGM in a density dependent fashion; de- 
creasing 7 increases the IGM temperature below the mean 
IGM density and decreases the temperature above mean 
density. If we assume H i to be in photoionisation equilib- 
rium, the neutral hydrogen density scales approximately as 
riHi oc rn'°'^r~,^A^~°'^'T~^' where Fm is the photoionisa- 
tion rate per H i atom (e.g. lRauchlll998h . Note that this is 
not the case for changes in To only, where the gas temper- 
ature (and hence also the residual H i fraction in the IGM) 
are altered independently of the gas density. 

These combined effects explain the changes seen in the 
PDF for different 7. As shown in Fig. (2] a lower value of 
7 decreases the amplitude of the PDF moderately over the 
flux range 0.1 < F < 0.9, an effect which is compensated for 
by a strong increase of the amplitude over the flux ranges at 
F — and F = 1. Absorption lines associated with moder- 
ately overdense regions become associated with colder, more 
neutral material and hence are narrower and have a higher 
opacity in the centre of the line. In contrast, underdense re- 
gions associated with transmission near the continuum level 
become hotter and therefore more highly ionised. Expressed 
in a different way, a larger 7 tends to erase the small scale 
power in the spectra by increasing the temperature of the 
gas associated with the most prominent absorption lines. 
Note that the effect of changing 7 is therefore partially de- 
generate with changing To and as within the observation- 
ally constrained ranges, although it has a significantly larger 
impact on the shape of the PDF. Fig. [5] suggests that a 
temperature-density relation which is close to isothermal is 
in better agreement with the K07 data. 

Finally, the effect of varying Tcs on the PDF is illus- 
trated in the right-most panels of Fig. (2] The value of Tcs 
measured in observed spectra exhibits rather large fluctua- 
tions between different lines of sight, and it is rather sensi- 
tive to the number of strong absorption systems (,Viel et al. 
l2004al ) and the contribution of metal absorption (K07). 
It is therefore reasonable to consider Toff as an indepen- 
dent free parameter within the measurement errors of the 
observations. The values displayed correspond to a mean 
flux range of ±0.020 around the fiducial values of (F) — 
[0.863,0.797,0.730] at (z) = [2.07,2.52,2.94]. The PDF is 
very sensitive to even small changes of Tea but in a quali- 
tatively different way to the other parameters. Altering Tcff 
changes the shape of the PDF by increasing the PDF am- 



plitude at one end and decreasing it at the other. This is in 
contrast to the effect of the other parameters, which when 
varied raise or lower the PDF at both ends simultaneously. 

3.3 Radiative transfer effects during reionisation 

Thus far we have assumed a tight, polytropic temperature- 
density relation for the low density IGM, T = ToA'^~^ , with 
1 < 7 < 1-6, which in the standard paradigm for the thermal 
state of the photo-ionised IGM is expected to apply to the 
majority of the optically thin gas responsible for the absorp- 
tion in the Lya forest. The relation arises from the balance 
between photo-heating and cooling by adiabatic expansion, 
with a slope which evolves from 7 ~ 1 following reionisa- 
tion to an asymptotic value of 7 — 1.6 (JHui fc Gnedinlll997l : 
iTheuns et aLlll998l : IValageas et al. ''2002Y Observational ev- 
idence for a boost in the IGM temperature and/or a flat- 
tening in the temperature-density relation a t 2: ~ 3 pro 



vides possible evidence for He 11 reionisation llSchave et al 
i^pO^; iRicotti et allbOOOl : iTheuns et a"Lll2002bl : IZaroubi et al 
I2OO6). As we saw in the last section, the observed flux PDF 
indeed appears to prefer a thermal state where the voids in 
the IGM are relatively hot, although this model is still not 
in perfect agreement with the data. It is therefore possible 
that a model with an inverted temperature-density relation 
with 7 < 1 may instead provide a better flt to the data. 

Is such a thermal state for the IGM plausible? During 
reionisation the assumption of a tight temperature-density 
relation for the low density IGM is likely to be a poor one. 
The filtering of ionising radiation by optically thick gas al- 
ters the spectral shape of the UV background, giving rise 
to additional heating effects which raise the gas tempera- 
ture and blur t he dependence of the IGM temperature on 
the gas de nsity ("Abel fc Haehnelt"l999': 'Bolton et al.l l2004l : 
iMaseUi fc F errara 2005; Tittley fc Mciksin 2007). The inho- 
mogeneous distribution of the ionising sources themselves 
will also co ntribute to spatial variations in the IGM ther- 
mal state llMiralda-Escude et all I2OO0I : iGleser et all l2005l : 
iFurlanetto fc OhI 120071 ). In addition, the cooling timescale 
for the low density gas responsible for the Lya forest is of 
order a Hubble time, meaning that changes to the thermal 
state of the IGM will persist fo r some time after reionisation 
(Mir alda-Escude fc Ree"3ll994h . 

Radiative transfer effects on the thermal state of the 
IGM may thus result in an inversion of the temperature- 
density relation, such that 7 < 1. As ionising photons prop- 
agate through optically thick, moderately overdense regions 
in the IGM, the spectral shape of the emission is consid- 
erably hardened as lower energy photons are preferentially 
absorbed. Any subsequent underdense gas is then subjected 
to ionising photons which have a substantially higher mean 
energy, boosting the temperature of the underdense region 
as it is photoionised and heated. One may envisage such a 
situation if ionising sources are initially embedded in over- 
dense, optically thick regions. State-of-the-art cosmological 
radiative transfer simulations indicate that this "inside-out" 
topology may be appropriate for modera te over and under- 
densities, at least during H i reionisation Ijlliev et al. 20061j). 
Such a radiative t ransfer induc e d tem perature inversion ef- 
fect was noted bv lBolton et al.l (|2004l ). although this result 
was a specific case for a single density distribution. More 
generally, radiative transfer effects during inhomogeneous 
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Figure 3. Contour plots of the volume weighted temperature-density plane at z = 2.07 in three of our hydrodynamical simulations. The 
number density of the data points increases by an order of magnitude with each contour level. The model in the left panel (20-256xs) has 
a temperature-density relation which exhibits a moderately increased scatter, while the central panel displays the simulation (20-256g3) 
with an inverted temperature-density relation. For comparison, the temperature-density plane in the fiducial model (20-256) is shown in 
the right hand panel. The dashed lines in each panel correspond to the temperature-density relations listed in Table 1. 



reionisation are expected to produce a more complex, mul- 
tiple valued relationship between temperature and density 
in the IGM. W hile we were in the final stages of preparing 
this manuscript iFurlanetto fc OhI (120071 ) presented a semi- 
analytical model aimed at describing the thermal state of 
the IGM during and after the inhomogeneous reionisation 
of helium. While probably not fully realistic due to a lack 
of a proper modeling of the filtering and spectral harden- 
ing effects of radiative transfer, the model confirms that 
an initially inverted and multi-valued temperature density- 
relation is indeed plausible. The exact character of these ef- 
fects will also depend on the types of sources which reionise 
the IGM l|Tittlev fc Meiksinll2007l ). 

Consequently, adopting a temperature-density relation 
with 7 < 1 may mimic some of the different heating effects 
expected due to inhomogeneous radiative transfer through 
an optically thick IGM. One may expect many under- 
dense regions to exhib it significa ntly elevated temp eratures 
following reionisation (I Abel fc Hachnclt 1999: Bolton et al.l 



12004 iGleser et al.ll2005l : iTittlev fc MeiksinI l2007l ) . Helium 
in the IGM at z ~ 3 exhibits patches which are still op- 
ticall y thick to photons with energies in excess of 4 Ry 
(e.g. Jjakobsen et al.|[l99i : iHeap et al.ll2000l : iReimers etahl 
I2OO5I ). Radiative transfer effects on the Lya flux PDF may 
therefore be important at the redshifts we consider here, 
especially if the tail end of He 11 reionisation o ccurs at 
3.2 (e.g.[Sp ngaila 1998; Schavo ct al. 2000; Thc uns et al 



2002bl. iBernardi et al.li 2003: SI mll et al.ll2004l : iBohon et al 
20061 : iFaucher-Giguere et al., ,20071 ') . Furthermore, using a 



semi-analytical model for the Lya forest based on den- 
sity distributions drawn from hydrodynamical simulations, 
LBeckor ct al. (2007) noted that an inverted temperature- 
density relation improved their fit to the flux PDF measured 
with an independent set of high resolution QSO spectra at 
1.7 <z< 5.8. 

The inhomogeneous spatial distribution of ionising 



sources is, however, also expected t o increase the scatte r 
in the temperature-density relation (JTheuns et al.ll2002d^ . 
We therefore compare two further different models to the 
K07 data, one with an inverted temperature-density relation 
(20-256g3) and another with a moderately increased scatter 
in the temperature-density relation (20-256xs). For illustra- 
tion, the volume weighted temperature-density planes for 
these two models ai z — 2.07 are displayed in Fig. [3] The 
flducial model is shown in the right-most panel. The dashed 
lines in each panel correspond to the temperature-density 
relations listed in Table 1. 

Fig. |4] displays the comparison between the flux PDFs 
constructed from these models to the observational data. 
Note that we have again matched t^s in the simulated spec- 
tra to the values measured by K07. The PDFs for the fidu- 
cial model (solid curves) and the model with increased scat- 
ter in the temperature-density relation (dashed curves) are 
very similar. At all redshifts the fiux PDFs exhibit differ- 
ences of only a few per cent. A moderate increase in the 
scatter of the temperature-density relation appears unlikely 
to substantially affect the shape of the PDF. However, the 
temperature-density relation is still relatively well defined 
in this model. There is also no good reason to assume a 
Gaussian dispersion. In reality the scatter may be rather 
more inhomogeneous. We therefore should not rule out the 
possibility that more complex, multiple valued relationship 
between temperature and density fits the K07 data better. 
We now turn to the model with an inverted temperature- 
density relation (dotted curves). In this case the modifica- 
tion of the thermal state of the IGM makes a substantial 
difference to the shape of the PDF, as expected from the be- 
haviour seen in Fig. [2] The inverted temperature-density re- 
lation brings the simulated and observed flux PDFs into ex- 
cellent agreement. The previously most problematic points 
at the extreme ends of the PDF at flux levels close to zero 
and unity now show an excellent match between the simu- 
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Figure 4. Comparison of the observational data of K07 at (2) = [2.07, 2.52, 2.94] (open diamonds with error bars) to simulations with 
an inverted temperature-density relation (20-256g3, dotted curve) and a model with increased scatter (20-256xs, dashed curve). The solid 
curve in each upper panel shows the PDF derived from the fiducial model (20-256), while the lower panels display the ratio of the other 
PDFs to the fiducial model. The shaded grey region in each panel corresponds to the portion of the PDF which is least affected by noise 
and continuum uncertainties. 



lations and data. Note again, however, that the IGM may 
nevertheless have a more complex thermal state following 
He II reionisation, rather than actually having a well de- 
fined, inverted temperature-density relation. We now pro- 
ceed to discuss the quality of the fits of different models 
more quantitatively in the next section. 



4 DETERMINATION OF THE BEST FITTING 
MODELS TO THE K07 PDF 

4.1 A quantitative comparison of the simulated 
and observed data 

As should be clear from the previous section and analyses 
by other authors, there is considerable degeneracy in the 
dependence of the flux PDF on cosmological and astrophys- 
ical parameters. The quality of the observational data and 
the apparently good agreement with the flux PDF from our 
inverted temperature-density relation simulation motivated 
us, however, to perform a quantitative x^ analysis. For this 
we used the full covariance matrix derived from the obser- 
vational data of K07. 

As discussed by iMcDonald et al.l (|2000l ). using the co- 
variance matrix is important; the errors bars on the PDF 
are highly correlated. However, the covariance matrices de- 
rived by K07 (their equation 1) are rather noisy for this 
purpose, preventing relia ble inversio n . We therefore follow 
the method suggested bv lLidz et al.l (|2006l ). where the ob- 
served covariance matrix is regularis ed using the corre lation 
coefficients, rs(i,j) = C0VB{i, j) / y^ cova{i, i)coVs{j , j) com- 
puted for each set of synthetic spectra, such that cov^ {i,j) = 
'''si'i, j) \/ covdii, i)covd{j, j) ■ Here covd and cova are the co- 
variance matrices computed from the observed data and 
synthetic spectra, respectively. Note that this implicitly as- 



sumes the observed and simulated spectra have similar co- 
variance properties. 

We perform a simple x^ analysis on the K07 data for 
two different cases. In the first instance we adopt a conserva- 
tive approach and analyse the PDF only over the fiux range 
0.1 < F < 0.8, where the impact of noise properties and con- 
tinuum uncertainties on the K07 PD F data are small. This 
is similar to the approach adopted bv lDesiacgues fc Nusserl 
(2005). However, as demonstrated in the last section, the 
extreme ends of the PDF at low and high flux levels contain 
important information on the thermal state of the IGM. For 
high resolution, high S/N data like that of K07, the uncer- 
tainties are expected to be relatively small over the entire 
range of the flux PDF. The continuum flts are generally ac- 
curate to within 1 — 2 per cent, and the typical S/N ~ 50 
(see K07 for a detailed discussion regarding the effect of 
these parameters on the PDF). We therefore also perform 
the x^ analysis over the entire PDF, < F < 1. 

In our analysis we do not attempt a global minimisation 
of x^ over the parameters we have varied in the simulations. 
We instead just tabulate the minimum reduced x^ values 
for our grid of models with different To, as, and 7. We vary 
Toff over an appropriate range, 25 per cent either side of the 
central values from the fit to the redshift evolution of the 
metal cleaned Lya effective optical depth reported by K07, 
Toff = 0.0023 ± 0.0007(1 + 2)3e5±02i. This scatter reflects 
the variance in the data. Continuum fltting uncertainties, al- 
though small at the redshifts we consider here, are another 
source of error in constraints on this quantity. We shall dis- 
cuss this in more detail shortly. 

The results of the x^ analysis for the flux range 0.1 < 
F < 0.8 are summarised in Table 3. In addition to the fldu- 
cial model (20-256) we tabulate the results for five other 
simulations with varying To, 7 and as- The models which 
have a low as, high Tq, high 7 and increased scatter in the 
temperature-density relation are omitted from Table 3, as 
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Table 3. The minimum x^ values for the best fitting models to the K07 PDF data at {z) = [2.07,2.52,2.94], fitted in the range 
0.1 < F < 0.8 only. The columns list the best reduced x^ (for f = 15 — 4 = 11 degrees of freedom) and the corresponding best fitting 
value for the effective optical depth, r^t. This is expressed as a ratio with the mean value of T^ff in each observational rcdshift bin. For 
reference these are r^g = [0.147, 0.227, 0.315] at (z) = [2.07, 2.52, 2.94]. 



Model 



{z = 2.07) 



(z = 2.52) 



(z = 2.94) 



Tflt/Tcff 

{z = 2.07) 



Tflt/feff 

{z = 2.52) 



Tflt/Toff 

(z = 2.94) 



Notes 



20-256 


4.57 


2.84 


2.03 


0.840 


0.857 


0.817 


Fiducial model 


20-256tl 


3.61 


2.49 


1.86 


0.864 


0.871 


0.822 


Low To 


20-256s2 


4.05 


2.65 


1.65 


0.886 


0.851 


0.819 


(T8 = 0.95 


20-256g2 


2.34 


1.90 


1.34 


0.898 


0.927 


0.849 


Low 7 


20-256s2g2 


1.85 


1.81 


1.14 


0.888 


0.925 


0.923 


as = 0.95, Low 7 


20-256g3 


1.62 


1.44 


0.90 


0.935 


1.025 


1.013 


Inverted 7 



Table 4. The x^ analysis results for the best fitting models to the whole PDF, < -F < 1. As in Table 3, the columns list the best 
reduced x^ (for i/ = 21 — 4 = 17 degrees of freedom) and the corresponding best fitting value for the effective optical depth, rat, as a 
ratio of the mean observed value in each redshift bin. For reference these are T^g = [0.147,0.227,0.315] at (2) = [2.07,2.52,2.94]. 



Model 



xV^ 

{z = 2.07) 



xV^' 

(z = 2.52) 



(z = 2.94) 



Tfit/Tcfi 

[z = 2.07) 



Tflt/Teff 

(2 = 2.52) 



Tfit/Tcff 

(2 = 2.94) 



Notes 



20-256 


11.15 


15.28 


13.47 


0.781 


0.833 


0.823 


Fiducial model 


20-256tl 


10.02 


15.22 


14.14 


0.808 


0.840 


0.844 


Low To 


20-256s2 


6.90 


9.65 


9.62 


0.820 


0.850 


0.854 


as = 0.95 


20-256g2 


5.03 


8.14 


8.03 


0.884 


0.914 


0.902 


Low 7 


20-256s2g2 


3.07 


4.70 


5.32 


0.923 


0.930 


0.927 


as = 0.95, Low 7 


20-256g3 


3.06 


1.56 


1.82 


1.057 


1.029 


1.025 


Inverted 7 



all have reduced x^ values which are similar to or worse 
than the fiducial model. The inverted temperature-density 
relation model (20-256g3) provides the best fit to the K07 
data with xV^^ = [1-62, 1.44, 0.90] at (z) = [2.07, 2.52, 2.94] 
for i/ = 11 degrees of freedom, corresponding to probabilities 
of [9,15,54] per cent. These are therefore reasonable fits in 
a quantitative sense. We would, however, urge the reader to 
exercise some caution in interpreting the absolute x^ values; 
estimating the observational PDF errors and the covariance 
matrix is problematic. We may have underestimated the ob- 
servational errors, especially at low flux levels. Note that the 
non-linear high density regions responsible for the absorp- 
tion at low flux levels are also the most sensitive to galactic 
winds and other uncertain details of the numerical simula- 
tions. It is nevertheless gratifying that the best fitting values 
for Teff are all within 7 per cent of the mean values in each 
redshift bin, well within the variance of this quantity. The 
model with an almost isothermal temperature-density rela- 
tion and as = 0.95 (20-256s2g2) is also in marginal agree- 
ment with the data. Note that although some of the remain- 
ing models are in reasonable agreement with the observed 
PDF for the highest redshift bin (z) — 2.94, they are all in 
poor agreement if the full redshift range of the K07 data set 
is taken into account. 

In Table 4 we list the minimum reduced x^ values for a 
comparison of the simulated and observed spectra over the 
full range of the PDF, < F < 1. The inverted temperature- 
density relation again provides the best fit to the K07 data 
with x'^/i' = [3.06,1.56,1.82] at (z) = [2.07,2.52,2.94] for 
1/ = 17 degrees of freedom. The formal probabilities of the 
reduced x^ values are now rather low, with a very small 
probability at {z) — 2.07 and [7,2] per cent probability at 
{z) — [2.52,2.94]. However, as already discussed, estimat- 



ing the errors of the PDF is difficult. Rather than being 
evidence against this model this may suggest that we have 
somewhat underestimated the errors at the extreme ends of 
the PDF, especially at low fiux levels. The best fitting val- 
ues for Toff are again within 6 per cent of the mean values 
in each redshift bin. All the other models are very bad fits. 

In Fig. \5\ we show the best fits to the K07 fiux PDF 
over the range 0-1 < F < 0.8 (solid curves) and < -F < 1 
(dotted curves). The models are almost indistinguishable to 
the eye in the two highest redshift bins, although one may 
clearly perceive the difference between the different PDFs 
at {z) = 2.07. The poor fit to the K07 data at (z) = 2.07 for 
< F < 1 is mainly due to the infiuence of the data point 
at F = 0; the agreement between the the observed and sim- 
ulated PDF at F > 0.8 is very good. Considering the results 
presented in Fig. [T] a model with strong galactic winds may 
also alleviate some of this discrepancy. Even though we have 
attempted to accurately model the noise properties of the 
observed spectra, differences in the noise level in the spec- 
tra may also play a small role in this disagreement (K07). 
In any case, the improved agreement between the observa- 
tional data and the Lya forest simulations with the inverted 
temperature-density relation is impressive. 

Finally, we note that the best fitting Tofi values for 
the inverted temperature-density relation model corre- 
spond to a photoionisation per hydrogen atom of Fhi = 
[1.43,1.05,0.97] X 10-^2 s'^ for 0.1 < F < 0.8 and 
Phi = [1.14,1.05,0.95] x 10"^^ s~^ for < F < 1 at 
{z) = [2.07,2.52,2.94]. These values are in excellent agree- 
ment with the results of lBolton et al.l (|2005r ) which are based 
on independent m easurements of the Lya forest opacity 
(|Schave et al.ll2003l ). 
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Figure 5. The best fit to the PDF data of K07 (open diamonds with error bars) for the inverted temperature-density relation model. The 
solid curves show the best fits over the flux range 0.1 < F < 0.8, corresponding to the data in the shaded portion of the plots, while the 
dotted curves show the best x^ fits to the entire PDF. The best fitting value of the effective optical depth for the models are also displayed 
in each panel. 



4.2 The effect of continuum uncertainties on the 
X^ analysis 

We now consider how important the continuum placement 
on the observed spectra is for the quality of the fit of the 
flux PDF at < F < 1. K07 discussed the effect of varying 
the continuum on the observational data in detail. For their 
typical estimated continuum uncertainty of ~ 1 — 2 per cent, 
the changes in the PD F are sm all and are most evid ent at 
F > 0.8. Tvt ler et al.1 (|2004) and lKirkman et all ifjobs') also 
estimated continuum uncertainties to be of order 2 per cent 
based on an analysis of a lar ge set of moderate res o lution 
Lya forest spectra at z < 3.2. lFaucher-Giguere et al.l (|2007|) 
estimate the continuum placement to be somewhat larger 
with a systematic bias towards an underestimate of the con- 
tinuum level by as much 4 per cent at z = 3 and 1 per cent 
at z = 2 in a recent measurement of the Lya effective op- 
tical depth using high-resolution data. Note, however, that 
while a systematic misplacement of the continuum appears 
likely with decreasing mean flux level (and thus increasing 
redshift) it will also depend strongly on the shape of the 
Lya flux PDF at high flux levels and thus on the tempera- 
ture of the gas in low-density regions of the IGM. 

We have checked for the presence of such a potential 
bias by comparing all our Lya forest models to the K07 data 
with synthetic spectra for which we have raised or lowered 
the continuum by 1.5 and 3 per cent. Raising the continuum 
on the synthetic spectra signiflcantly worsens the agreement 
between the K07 data and all the simulations. Lowering the 
continuum on the synthetic spectra, mimicking a contin- 
uum which has been placed too low on the observational 
data, also strongly increases the x^ values, although the 
models with either an inverted or isothermal temperature- 
density relation are still slightly preferred in some redshift 
bins. Misplacing the continuum changes the flux PDF in a 
complex way. The fact that the agreement between the ob- 
served data and all models worsens significantly on varying 
the continuum suggests that continuum fitting errors for the 
high-quality data of K07 are under control at a level of 1- 
2 per cent, consistent with the estimates of the continuum 



uncertainties by K07. Altering the continuum placement on 
the synthetic spectra does not alter our conclusions regard- 
ing the best fitting model to the K07 data. 



4.3 Comparison to other studies 

Several other detailed studies of the Lya fiux PDF are pre- 
sented within the existing literature, and we now discuss 
som e of these results i n com parison to this work. 

iMcDonald et al.1 (|2000l ') (hereafter MOO) presented a 
measurement of the flux PDF at (z) = [2.41, 3.00, 3.89] from 
a sample of eight high resolution QSO spectra. A detailed 
comparison between the widely used MOO measurement and 
the measurement used in this work is given in K07. Here 
we comment on the interpretation of the MOO data when 
compared to hydrodynamical Lya forest simulations. MOO 
compare their observational data to the LIO Eule ri an hy - 
drodynamical simulation of iMiralda-Escude et al.l ([199g). 
The simulation parameters were ilm = 0.4, Qa = 0.6, 
f^b/i^ = 0.015, h = 0.65, as = 0.79, n = 0.95, To = 16 000 K 
and 7 ~ 1.3 at 2 = 3, with a box size of 10h~^ comoving Mpc 
and a 288"^ grid. Performing a x^ analysis with the covari- 
ance matrix, MOO obtain x^/'^ = [4-4, 1.7, 1.1] for 19 degrees 
of freedom at (2) = [2.41,3.00,3.89]. Note, however, that 
they do not minimise over TeS as we do. Instead, MOO rescale 
the synthetic spectra to match the measured effective optical 
depth of the observational data, t^s = [0.201,0.380,0.744]. 
The fit is performed over the range < F < 1, and they 
attribute the disagreement at (z) = 2.41 as being due to the 
effect of continuum uncertainties on the PDF. 

As discussed in the last section this is in contrast with 
our findings: using updated cosmological parameters we 
find the K07 PDF, which is systematically lower at 0.2 < 
F < 0.7 than the MOO data, strongly favours an inverted 
temperature-density relation when fitting over < F < 1. 
It is not obvious if this difference is due to the different 
shape of the measured PDFs, differences between the model 
parameters or differences in the numerical simulations. To 
investigate this further we have attempted to fit the fiux 
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PDF of our simulations ai z — 2.94 to the MOO measure- 
ment of the flux PDF at (z) = 3. The small difference in 
redshift should not be important here. When minimising 
over TeB we find that all of the models listed in Table 3 are 
in good agreement with the MOO data for 0.1 < F < 0.8 (all 
have X^/i^ < 1.04). The inverted temperature-density rela- 
tion model still provides the best reduced x^ ix^ /'^ ~ 0-41 
for Tcff — 0.411) perhaps suggesting that the errors in this 
flux range have been overestimated by MOO. When perform- 
ing the x^ analysis over the entire MOO PDF, we were at 
first not able to obtain an acceptable fit to the MOO data. 
Note, however, that MOO also adjusted the normalisation of 
their synthetic spectra to take into account possible contin- 
uum errors of the observational data. They achieved this by 
identifying the maximum transmitted fiux in each synthetic 
line-of-sight, F^, as the continuum flux and then proceeded 
to divide the flux in all other pixels in the line-of-sight by 
Fc 

We therefore again raise or lower the continuum in our 
synthetic spectra by 1.5 per cent, this time for comparison 
to the MOO data at (z) = 3. As for the K07 data, raising 
the continuum on the synthetic spectra signiflcantly wors- 
ens the agreement between the MOO data and all our models. 
Lowering the continuum in this instance, however, improves 
the agreement between the MOO data and simulations re- 
markably. Assuming our synthetic spectra provide an accu- 
rate representation of the Lya forest, this suggests that the 
continuum in the absorption spectra on which the MOO mea- 
surement is based may indeed be placed somewhat low. If we 
lower the continuum level by 1.5 per cent, we flnd the model 
with a temperature-density relation close to isothermal and 
a larger as — 0.95 (20-256s2g2) provides the best fit to the 
MOO data, with x^/'^ = 1-08 for 8-n effective optical depth 
TeS = 0.385, which is pleasingly close to the MOO value. The 
inverted temperature-density relation model (20-256g3) pro- 
vides the next best fit with x^/^ = 1-91 for Tefi = 0.418. This 
is in contrast to MOO, who achieved a reasonable fit to the 
data using different (now somewhat outdated) cosmologi- 
cal parameters and a standard temperature-density relation. 
The different conclusion we reach in this work is therefore 
probably largely due to the updated cosmological param- 
eters and perhaps also numerical differences in the simula- 
tions. In summary, we conclude that the MOO PDF measure- 
ments are consistent with our numerical simulations if the 
continuum is lowered by 1.5 per cent. The MOO measurement 
also appears to favour a flattened or inverted temperature- 
density relation at z ~ 3, consistent with our results based 
on the independent K07 data set. 

A more recent joint analysis of the best fltting model 
to the MOO PDF and flu x powe r spectrum was undertaken 
by iDesiacaues fc Nusseij (|200i) at z = 3. They used N- 
body simulations with 256'^ dark matter particles and a box 
size of 25h^^ comoving Mpc. Synthetic Lya forest spectra 
were constructed by smoothing the dark matter density and 
velocity fields with a Gaussian filter over a characteristic 
scale, which was left as an extra free parameter in their 
analysis. Unfortunately, this extra free parameter compli- 
cates our comparison to this work as it is degenerate with 
other relevant quant i ties, most notably 7. Consequently, 
IDesiacaues fc Nusserl (|2005l ) d o not advocate a value f or this 
quantity. On the other hand, IDesiacaues fc Nusserl (120051 ) 
also found that their joint analysis of the MOO flux power 



spectrum and PDF (fitted over the range 0.05 < F < 0.75) 
may favour a lower valu e of as than previous analyses of the 
power spectrum alone (IViel et al.l l2004al : iMcDonal d et al.l 
12003' ). although this result should again be degenerate with 
the thermal state of the IGM. We instead find that the 
K07 PDF measurements on their own favour value of as 
towards the upper end of current observational constraints 
if t he temperat u re-den sity relation is close to isothermal. 

iLidz et al.l l|2006l ) compare synthetic spectra drawn 
from pseudo-hydrodynamical simulations to the MOO ob- 
servational data, but instead use a different estimator for 
the flux PDF. This estimator is derived from data which is 
smoothed on two different scales, thereby having the advan- 
tage of rendering the observational data insensitive to the 
normalisation and shape of the continuum, but also to struc- 
ture on small scales in the spectra. This latter effect may be 
advantageous when using simulations which are not fully re- 
solved, but it may also erase important information in the 
observational data , especially with regard to the thermal 
state of the IGM. iLidz et al.l l|2006l ) find their simulations 
favour a temperature-density relation with 7 > 1.32 {2a) 
a,t z — 2.72. However, given th e approxim a te tre atment of 
smoothing due to gas pressure iLidz et al.l (|2006l ) adopt in 
their simulations, as well as their different observational es- 
timator for the PDF, the significance of this constraint on 
comparison to our results is difficult to judge. 

Lastly, and as discussed previously, iBecker et al.l (|2007l ) 
found that an inverted temperature-density relation may be 
required to adequately fit the flux PDFs measured from an 
independent set of high resolution QSO spectra over a wide 
redshift range, 1.7 < 2 < 5.8. Our results ar e also consis- 
tent w ith this possibility. Note, however, that iBecker et al.l 
|2007J) use a semi-analytical model for the Lya forest based 
on the IGM density distributions drawn from a hydrody- 
namical simulation. This model does not self-consistently 
model Jeans smoothing effects, and a direct comparison to 
their results is again not possible. 



5 IMPLICATIONS FOR THE Lya FOREST 
FLUX POWER SPECTRUM 

The Lya fiux PDF is of course not the only diagnostic of the 
thermal state of the IGM. There should be good prospects 
to break some of the degeneracies discussed here and fur- 
ther tighten constraints by a combined analysis with other 
Lya flux statistics. Assumptions about the IGM thermal 
state are of particular relevance for measurements of the 
matter power spectrum from the Lya flux power spectrum. 
Such analyses have so far not allowed for the possibility of an 
inverted temperature-density relation and have marginalised 
over what was believed to be the plausible range for the 
slope, 1 < 7 < 1.6. Due to the temperature dependence of 
the recombination coefficient a smaller 7 should translate 
into a steeper power law relation between neutral hydrogen 
density and matter density. For measurements of the matter 
power spectrum a smaller assumed 7 should therefore result 
in a smaller amplitude for the matter power spectrum. In 
this section we therefore briefly discuss the implications of 
an inverted temperature-density relation for the Lya forest 
flux power spectrum. Note this is a qualitative comparison 
only; a detailed joint analysis of the Lya fiux PDF and flux 
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density relation therefore mimics some of the effects a larger 
normalisation for the matter power spectrum heis on the 
Lya forest flux power spectrum. 

This has potentially important implications for infer- 
ring the matter power spectrum from the Lya forest. The 
measurements of the matter power spectrum from Lya for- 
est data are affected by non-linearities as well as by complex 
degeneracies between the thermal state of the IGM, the den- 
sity distribution and the assumed effective opti cal depth. 
Howeve r, if we naively extrapolate the results of IViel et al.l 
((2004b|) to smaller values of 7 we would expect a 5-10 per 
cent percent lower amplitude ag for the slope of 7 ~ 0.5 
favoured by the analysis of the Lya flux PDF presented 
here. A 5-10 per cent lower amplitude for the inferred matter 
power spectrum would largely alleviate the strain between 
the measureme nts from Lya forest data and the WMAP 3 
year results (Viel&JHaehiiclt 2006). Unfortunately it would 
also substantially weaken the strong upper limits on the neu- 
trino mass from the joint analysis of Lya for est and CMB 
data (jSeliak et al.ll2006al : iGratton et al.ll2007l ). 



Figure 6. The ratio of the Lya flux power spectrum computed 
from synthetic spectra drawn from the models indicated on the 
panel at 2 = 2.94 to the flux power spectrum of the fiducial model 
(20-256, ITS = 0.85, 7 = 1.27). The shaded region corresponds to 
the range of wavenumbers used by Viel et al. (2004b) to infer 
the amplitude and shape of the matter power spectrum from the 
Lya flux power spectrum, 0.003 < fc(s km~^) < 0.03. An in- 
verted temperature-density relation increases the power within 
this wavenumber interval, mimicking the efl'ect of a larger ampli- 
tude (Tg of the underlying matter power spectrum. 



power spectrum with a larger suite of numerical simulations 
with inverted temperature-density relations will be required 
to examine this more accurately, which is beyond the scope 
of this paper. 

Fig. |6] illustrates the effect of an inverted temperature- 
density relation on the flux power spectrum at z = 2.94. 
The solid curve shows the ratio of the flux power spec- 
trum, PF(fc), computed from the model with the inverted 
temperature-density relation to the fiducial model power 
spectrum, PF,2o-25eik). The shaded region cor r espond s to 
the range of wavenumbers used by I Viel et al.l (|2004bj) to 
infer the matter power spectrum from the Lya forest flux 
power spectrum. The flux power spectrum of the simula- 
tions with the inverted temperature-density relation clearly 
exhibits more power over this range on comparison to the 
flducial model, and this increases noticeably towards larger 
wavenumbers. As discussed in section 13.21 this is because 
a larger 7 increases the temperature and hence the ionised 
fraction of the gas associated with the most prominent ab- 
sorption features, thereby suppressing predominantly the 
small scale power. The dotted curve shows the same compar- 
ison for the isothermal temperature-density relation model. 
The effect is again to increase power at all scales, although 
to a somewhat lesser extent. The dashed and dot-dashed 
curves also show the comparison for two models with differ- 
ent matter power spectrum normalisations, as — 0.75 and 
as ~ 0.95, respectively. As expected, a lower as suppresses 
power and a larger as increases power over the relevant 
wavenumber range. An inverted or isothermal temperature- 



6 CONCLUSIONS 

We have compared the improved measurements of the 
Lya flux PDF at 1.7 < 2 < 3.2 by K07 obtained from a fuU 
Voigt profile analysis of 18 high resolution, high signal to 
noise UVES/VLT spectra to a large set of hydrodynamical 
simulations of the IGM. Our main findings are as follows: 

• The shape of the Lya flux PDF depends only weakly 
on the fluctuation amplitude of the matter power spectrum 
characterised by as and on the normalisation of the 
temperature-density relation for the low density IGM , 
T = ToA^-i (.Hui fc Gne dh] 1 19971 : IValageas et al.l |2002| ). 
A moderately increased scatter in the temperature-density 
relation, which may result from inhomogeneous radiative 
transfer effects during the epoch of He 11 reionisation, also 
has little effect on the shape of the flux PDF. The shape 
of the Lya flux PDF is mainly determined by the effective 
optical depth Tcff and by the slope of the temperature- 
density relation for the low density IGM. The extreme ends 
of the PDF at flux levels of zero and unity are particularly 
sensitive to the slope of the temperature-density relation. 

• A comparison of the full range of the Lya flux 
PDF to hydrodynamical simulations favours an inverted 
temperature-density relation. The flux PDFs of our syn- 
thetic spectra with an inverted temperature-density relation, 
such that 7 = [0.67,0.54,0.44] at (z) = [2.07,2.52,2.94], 
provide an excellent fit to the observed metal-cleaned 
flux PDF measured by K07. This suggests that the voids 
in the IGM are significantly hotter than previously assumed. 

• If we restrict our analysis of the PDF to flux levels 
0.1 < F < 0.8, which are least affected by continuum 
fitting, noise, galactic winds and the uncertainties of how 
to simulate the relevant physics in dense regions of the 
IGM, the comparison of observed and synthetic spectra still 
favours a model where the temperature-density relation is 
either nearly isothermal or inverted. 
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• Raising or lowering the assumed continuum level of 
our synthetic spectra in order to mimic the possibility of 
a systematic over or underestimate of the continuum level 
in the observed spectra rapidly worsens the quality of the 
fit to the K07 PDF for offsets larger than 1.5 per cent. 
This suggests that continuum fitting errors for the high 
quality data of K07 are under control at a level of 1-2 
per cent, consistent with the estimates of the continuum 
uncertainties by K07. 

• A comparison of our simulations with the full range 
of the Lya flux PDF measured by MOO at (z) = 3 from 
a smaller observational sample gives consistent results 
if the continuum on the synthetic spectra is lowered by 
1.5 per cent. Using the currently preferred cosmological 
parameters, the MOO data also favour a nearly isothermal 
or inverted temperature-density relation. 

• An inverted temperature-density relation has important 
implications when inferring the underlying matter power 
spectrum from the Lya forest flux distribution. The effect of 
an inverted temperature-density relation on the flux power 
spectrum is similar to the effect of an larger normalisa- 
tion, (Tg, of the matter power spectrum. This could allevi- 
ate some of the existing strain between the measurements 
of erg from Lya f orest data and the WMAP 3 year results 
l|Viel fc Haehneh|[2006; ). 

The preference of the data for an inverted temperature- 
density relation would, if consolidated, be a rather profound 
result. It would suggest that the widely accepted paradigm 
of a tight IGM temperature-density relation at low densities, 
evolving from being almost isothermal following hydrogen 
reionisation to a r elation w here the temperature rises with 
increasing densitv (|Hui fc G ncdin 1997; Thouns ot al. 19981: 
IValageas et al.ll2002r ). is too simplistic. 

Perhaps the most plausible explanation for an inverted 
temperature density relation is He ii reionisation. There is 
evidence to suggest this may end around z ~ 3 (ISongaila 



19981 : ISchave et all l2000l : iTheuns et al.l r2002al: IShull et al 



2004 : iBolton et aiTl2006l : iFaucher-Giguere et al.lT2007l ). and 
thus lies at the upper end of the redshift range considered 
in this work. Changes in the photoheating rate which oc- 
cur as ionising photons propagate through optically thick 
regions can signiflcantly raise the temperature of the IGM 
l|Abel &: Haehnehi [l999|). First attempts at modeling in- 
homogeneous He II reionisation and/or the associated ra- 
diative transfer effects also suggest a more complex and 
possibly multi-value d and inverted temperature-dens i ty re- 
lation is plausible llBolton et al] l2004l: [cTeser et all [20051 : 



iTittlev fc MeiksinI [20071 : iFurlanetto fc Ohll2007l ). Therefore, 
while a well-defined, inverted temperature-density relation 
may not faithfully represent the true thermal state of the 
IGM, a model with 7 < 1 should mimic some of the ex- 
pected inhomogeneous radiative transfer induced effects on 
the IGM temperature following He 11 reionisation. Within 
this framework, this model provides a good fit to the K07 
observational data and suggests that the thermal state of the 
IGM may be significantly more complex at 2 ~ 3 than pre- 
viously thought. Further investigation of the thermal state 
of the IGM around z ~ 3 is therefore highly desirable. Note, 
however, that studies using wavelets have failed to detect 



the expec ted fluctuations in the IGM temperatu re at these 
redshifts (JTheuns et al.ll2002dl : [Zaldarriagall2002l '). 

While radiative transfer effects and the inhomogene- 
ity of He 11 reionisation provide a plausible explanation 
for an increased temperature in underdense regions of the 
IGM, it will require full numerical radiative transfer sim- 
ulations to conflrm the validity of this picture. Helium 
reionisation is driven by hard photons emitted by rare, 
luminous and most likely short-lived sources (QSOs) un- 
der conditions where photon mean free paths are signifl- 
cantly longer and recombination time scales shorter (if com- 
pared to the Hubble time) than for hydrogen reionisation at 
higher redshifts. Current 3D cosmological radiative trans- 
fer codes which incorporate multiple sources are designed 
primarily to investigate the topology of H i reionisation, 
and are less reliable w hen mode ling the resulting ther- 
mal state of the IGM l|lliev et al . 2006a). Most detailed 
studies of radiative transfer effects on the IGM tem- 
perature therefore employ either reduced dim ensionality 
(JBolton et al.l |2004| : iBolton fc Haehneltl 12007311 or neglect 
the modeling of multiple sources i Maselli fc Ferraral l2005l : 
Tittloy & Moiksin 2007). Furthermore, the physical scales 
involved during He 11 reionisation are much larger than 
those during H i reionisation at z > 6, placing even greater 
de mands on the required computationa l resources (although 
see Soka sian et al.ll2002l : IPaschoj 120071 ). Radiative transfer 
simulations which accurately track the temperature evolu- 
tion are thus very challenging. Developments in this area 
will be key to interpreting the increasing amount of high 
quality Lya forest data. 

Finally, we must also not lose sight of the possibil- 
ity that other pr ocesses may contr i bute to the heating 



of th e IGM (e.g. iNath et alj [l999l: iMadau fc Efstathioul 



19991: llnoue fc Kamaval l2003l : iRicotti fc Ostriked |2004| : 



Samuietal. l l 200,5h. O ;; r^^t. may instead indicate^^t 



an alternative modification to the current models for the 
Lya forest and the thermal state of the IGM is required. 
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